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Abstract 

A mean-field approach (filtering out subgrid scales) is applied to the Boltzmann 
equation in order to derive a subgrid turbulence model based on kinetic theory. 
It is demonstrated that the only Smagorinsky type model which survives in the 
hydrodynamic limit on the viscosity time scale is the so-called tensor-diffusivity 
model. Scaling of the filter-width with Reynolds number and Knudsen number is 
established. This sets the first rigorous step in deriving turbulence models from 
kinetic theory. 



1 Introduction 



The application of a filtering procedure to equations of hydrodynamics (Navier- 
Stokes equations) in order to construct a subgrid model is often used for the 
turbulence modeling [1]. The aim of such models is to take into account the 
effects of subgrid scales as an extra stress term in the hydrodynamic equations 
for the resolved scale fields. Further, the subgrid scale terms should be repre- 
sentable in terms of the resolved fields. This procedure, like any other attempt 
to coarse-grain the Navier-Stokes equations, runs into the closure problem due 
to the nonlinearity of the equation and due to the absence of scale separation. 
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On the other hand, in statistical physics, good schemes to obtain closure ap- 
proximations are known for nonlinear evolution equations (with a well-defined 
separation of scales). Unfortunately, attempts to borrow such schemes fail for 
the Navier-Stokes equations. The fundamental reason for this failure of the 
coarse-graining procedures on the Navier-Stokes equations is the absence of 
scale separation. Further, the length over which the equation is coarse-grained 
(the filter width in the present case) is completely arbitrary (and in practice 
dictated by the available computational resources), and cannot be justified a 
priori on physical grounds. 

In this paper, we show that a coarse-grained description of hydrodynamics 
using the microscopic theories is possible. Specifically, we apply the standard 
filtering procedure (isotropic Gaussian filter) not on the Navier-Stokes equa- 
tions but on the Boltzmann kinetic equation. We recall that the Navier-Stokes 
equations are a well defined limit of the Boltzmann equation (the hydrody- 
namic limit), whereas the filtering operation and going to the hydrodynamic 
limit are two distinct operations which do not commute, because kinetic fluc- 
tuations generally do not annihilate upon filtering. By doing so, we obtain the 
following results: 

• Smallness parameter: The smallness parameter of the present theory is the 
usual kinetic-theory Knudsen number Kn, 



where Ma is the Mach number and Re is the Reynolds number, v is the kine- 
matic viscosity, c s is the sound speed and L is the characteristic macroscopic 
length. Smallness of Kn rules emergence of both, the usual viscosity terms, 
and the subgrid contributions, on the viscosity time scale of the filtered 
Boltzmann equation (that is, in the first-order Chapman-Enskog solution 
to the filtered Boltzmann equation). 
• Scaling: In the coarse-grained representation obtained by filtering, the filter- 
width A (for the Gaussian filter, A 2 is proportional to the covariance) is 
the smallest length-scale to be resolved. The requirement that contributions 
from the subgrid scales appear in the kinetic picture at the time scale of 
molecular relaxation time (viscosity time scale) sets the scaling of A with 
the Knudsen number as follows: 



where k is a nonuniversal constant which scales neither with L, nor with 
Kn. For the sake of simplicity, we set k — 1 in all the further computations. 
Equations (1) and (2) imply that the filter-width scales with the Reynolds 
number as follows: 
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(2) 



A ~ Re" 1 / 2 . 



(3) 
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While the Kolmogorov length, l K , scales as Ik ~ Re 3 ^ 4 , we have 

^ ~ Re 1 / 4 . (4) 
Ik 

Thus, the filtering scale is larger than the Kolmogorov scale when Re is 
large enough. 

Subgrid model: With the above smallness parameter (1), and the scaling (2), 
we rigorously derive the following subgrid pressure tensor P^p, in addition 
to the usual (advection and viscosity) terms in the momentum equation: 
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Here p is the filtered density, and summation convention in spatial compo- 
nents is adopted. For any function X, X denotes the filtered value of X. 
Furthermore, the filtered rate of the strain tensor S a p and the filtered rate 
of the rotation tensor Q a p depends only the large scale velocity, 

l _ (6) 

5a/3 = - {O a Up + dpU a } . 

The derived subgrid model belongs to the class of Smagorinsky models 
[4], and the tensorial structure of the subgrid pressure tensor (5) corre- 
sponds to the so-called tensor- diffusivity subgrid model (TDSG) introduced 
by Leonard [5], and which became popular after the work of Clark et al 
[6]. Here, it is interesting to recall that in the class of existing Smagorin- 
sky models the TDSG is one of only a few models in which the sub-grid 
scale stress tensor remains frame-invariant under arbitrary time-dependent 
rotations of reference frame [7]. Furthermore, the TDSG model belongs to a 
subclass of Smagorinsky models which take into account the backscattering 
of energy from the small scale to the large scales [1]. Beginning with the 
seminal work of Kraichnan [8], importance of the backscattering of energy 
in turbulence modeling is commonly recognized. 

Uniqueness: The result (5) requires only isotropy of the filter but otherwise is 
independent of the particular functional form of the filter. There are no other 
subgrid models different from (5) which can be derived from kinetic theory 
by the one-step filtering procedure. In other words, higher-order spatial 
derivatives are not neglected in an uncontrolled fashion, rather, they are of 
the order Kn 2 , and thus do not show up on the viscosity time scale. 
Nonarbitrary filter-width: Unlike the phenomenological TDSGM where the 
prefactor in Eq. (5) remains an unspecified "A 2 ", kinetic theory suggests 
that the filter-width cannot be set at will, rather, it must respect the phys- 
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ical parameterization (specific values of Re, Kn etc) of a given setup when 
the subgrid model is used for numerical simulation. Recent findings that 
simulations of the TDSG model become unstable for large A [9] is in qual- 
itative agreement with the present result that the filter-width A cannot be 
made arbitrary large. 

The structure of the paper is as follows: In section 2 we set up the kinetic 
theory for the subsequent coarse-graining and derivation of the subgrid model. 
It is important to stress that the only requirement on the choice of the kinetic 
equation in the present context is that it gives the Navier-Stokes hydrodynamic 
equations in the appropriate fully resolved limit. For that reason we choose to 
work with a recently introduced minimal kinetic model [2,3] which is sufficient 
for our present purpose. Filtered kinetic equation is obtained in section 3. In 
section 4 we derive the subgrid model (5) using the Chapman-Enskog method 
[10] for the filtered kinetic model. This derivation uniquely defines the scaling 
(3) from the requirement that the subgrid terms appear on the viscous time 
scale. Finally, results are summarized, and some directions of future research 
are discussed in section 6. 



2 Kinetic theory 

For the present discussion, the particular choice of the kinetic model is unim- 
portant as long as the hydrodynamic limit of the kinetic theory is the usual 
Navier-Stokes equations at least up to the order 0(Ma 3 ). We demonstrate the 
whole procedure in detail for a recently introduced minimal discrete-velocity 
kinetic model [2,3]. As the final result is just the same in two and three di- 
mensional cases, for the sake of simplicity we chose to demonstrate the whole 
procedure using the two-dimensional model (D — 2). The kinetic equation is, 



The local equilibrium / 4 eq is the conditional minimizer of the the entropy func- 
tion H: 




(7) 



C x = {0,1,0,-1,0,1,-1,-1,1}, 
C y = {0,0,1, 0,-1, 1,1, -1,-1}. 



(9) 
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under the constraint of local conservation laws: 



EHi, Q}. = {p, P u}. (ii) 

i=i 



The weights in the equation (10) are: 

^ = ^{16,4,4,4,4,1,1,1,1}. (12) 
The explicit expression for /f q reads: 

r = m n (2 - n/TTw) ( 2 "-+/rW ) c - (13) 

Q=l ^ / \ 1 M Q / 



Below, it will prove convenient to work in the moment representation rather 
than in the population representation. Let us choose the following orthogonal 
set of basis vectors in the 9-dimensional phase space of the kinetic equation 
(7): 

^ = {1,1,1,1,1,1,1,1,1}, 

V> 2 = {o,i, 0,-1, o,i, -1,-1,1}, 

Vs = {0,0, 1,0, -1,1, 1,-1,-1}, 
V> 4 = {0,0, 0,0, 0,1, -1,1,-1}, 

ij 5 = {0,1, -1,1, -1,0, 0,0,0}, (14) 
^6 = {0,-2, 0,2, 0,1, -1,-1,1}, 
^ 7 = {0, 0,-2, 0,2, 1,1, -1,-1}, 
^ 8 = {4, -5, -5, -5, -5, 4, 4, 4, 4}, 
ij 9 = {4, 0,0, 0,0, -1,-1, -1,-1}. 

The orthogonality of the chosen basis is in the sense of the usual Euclidean 
scalar product, i.e., 

9 

J2^iktpkj = diSij, (15) 
k=i 

where d{ are some constants needed for the normalization (the basis vectors 
are orthogonal but not orthonormal). We define new variables Mi, i — 1, ... ,9 
as: 

U X '••<,//• (16) 

3=1 

where ipij denotes jth component of the 9-dimensional vector ipi. Basic hy- 
drodynamic fields are M\ = p, Mi = pu x , and M3 = pu y . The remaining six 
moments are related to higher order moments of the distribution (the pressure 
tensor P a p = fiCi a Cip and the third order moment Q a/ 3^ = Y fiCiaCipC^ 
and so on), as: M 4 = P xy , 
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M 6 = 3j2f i C? y C ix -2M 2 , 
M 7 = 3 £ WlC iy - 2M 3 



(17) 



The explicit form of the stress tensor in term of the new set of variables is: 



P xy — M 4 , 

Pxx = \Mi + ^M 5 + ^M 8 - ^Mg, (18) 
Pyy = \Mi ~ \m 5 + ^M 8 - ±M 9 . 



The time evolution equations for the set of moments are: 



d t M x + d x M 2 + d y M 3 = 0, 
d t M 2 + d x Qm x + ^M 5 + ^M 8 - Im 9 ) + <9,M 4 = 0, 

d t M 3 + d x M A + d y Qm, - l -M 5 + ^M 8 - ^M 9 ) = 



d t M 4 + ^ (2M 3 + M 7 ) + \d y (2M 2 + M 6 ) = - (M 4 cq (M 4 , M 2 , M 3 ) - M 4 ) , 
dot 

ftM 5 + ^ (M 2 - M 6 ) + ^ (M 7 - M 3 ) = -j- (M 5 cq (M 4 , M 2 , M 3 ) - il^g) 

d t M 6 - \d x (5M 5 - M 8 + M 9 ) + <9,M 4 = - (M 6 eq (M 4 , M 2 , M 3 ) - M 6 ) , 
o r 

d t M 7 + d x M A + \d y (5M 5 + M 8 - M 9 ) = - (M 7 eq (M 4 , M 2 , M 3 ) - M 7 ) , 
o r 

d t M s + d x (M 2 + 3M 6 ) + 0„ (M 3 + 3M 7 ) = - (M 8 cq (M 4 , M 2 , M 3 ) - M 8 ) , 

r 

d t M 9 - (2M 2 + M 6 ) - ^ (2M 3 + M 7 ) = - (M 9 eq (M 4 , M 2 , M 3 ) - M 9 ) . 



The expression for the local equilibrium moments M 4 eq , % — 4, . . . , 9 in terms 
of the basic variables M 1; M 2 , and M 3 to the order w 2 is: 
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M 6 eq (Mi,M 2 ,M 3 ) 
M^(M U M 2 ,M 3 ) 



Mf - Mf 



-M 2 , 
-M 3 , 



(20) 



m;\m 1 ,m 2 ,m 3 ) 



Mf + Mf 



MP(M 1 ,M 2 ,M 3 ) 




3 (M| + Mf 
Mi 



) 



The incompressible Navier-Stokes equations are the hydro dynamic limit of the 
system (19) and (20). 

In the next section, we shall remove small scales through a filtering procedure 
on the moment system (19), (20). A precise definition of the small scales is 
postponed until later sections. For the time being, let us assume that there 
exist a length-scale A, and we wish to look at the hydrodynamics at length- 
scale larger than A only. 



3 Filtered kinetic theory 

Coarse-grained versions of the Boltzmann equations have been discussed in the 
recent literature [11,12,13]. However, a systematic treatment is still lacking. 
In this section, we shall fill this gap. 

3. 1 Gaussian filter 

For any function X, the filtered function X is defined as: 



Function G is called the filter. In the sequel, we apply the filtering operation 
(21) on the moment system (19). We will need two relations. First, for any 
function X, 



This relation is sufficient to filter the propagation terms in the equation (19) 
due to linearity of propagation in the kinetic picture. The latter is a useful 




(21) 



d a X = d a X. 



(22) 
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property which is not shared by the hydrodynamic Navier-Stokes equations, 
where the nonlinearity and nonlocality both come into the same (uVn) term. 
Any isotropic filter, which satisfies the condition of commuting of the deriva- 
tives under the application of the filter (Eq. 22), will suffice for the present 
purpose. We choose a standard Gaussian filter [1] which has the property (22): 




G(r,A)=(-^) exp(-— ]. (23) 



Let us recall the isotropy properties of a Gaussian filter: 



/ G(r,AW = l, 
Jr d 

J G(r,A)rdr = 0, (24) 

r A 2 
Jrd G ^ A ^ r ' ar P dr = u 6a ?- 



Second, in order to filter the nonlinear terms (20) in the right hand side of 
moment equations (19), we will also need the following relation for three arbi- 
trary functions X, Y, Z which follow immediately from the isotropy property 
by second-order Taylor expansion: 



(cf) = i^r) + ill ^ x ^ d ^ - h {daZ) ( xd « Y + Yd « x + 2 -§~ daZ " 



+ 0(A 4 



(25) 



The effect of a Gaussian filter need not be truncated to any order at the 
present step. The higher-order terms lumped under 0(A 4 ) in equation (25) 
can be computed from elementary Gaussian integrals. As we shall see it soon, 
higher than second order terms disappear in the hydrodynamic limit once the 
scaling of the filter-width versus Knudsen number is appropriately chosen. 

In the next section, we shall filter the moment equations (19). 



3.2 Filtering the moment system 



Applying the filter (21) to the moment system (19), (20), using (22) and (25), 
and keeping terms up to the order u 2 , we obtain the following filtered moment 
system: 
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d t M 1 + d x M 2 + d y M 3 =0, 

0,JT, 0, (|mx + hf 5 + ^M 8 - l -M^j + d y M 4 = 0, 
'2^ 1^ 1^ 1- 



d t M 3 + a,M 4 + 9, (-M! - -M 5 + -M 8 - -M 9 ) =0, 

9 t M 4 + ^ (2M 3 + M 7 ) + ^ (2M 2 + M 6 ) =- (m^Mx, M 2 , M 3 ) - M 4 ) 



3 

+ T2Tm7^^ + (v 



d t M 5 + ^d x (M 2 - M 6 ) + (M 7 - M 3 ) = l - (M 5 eq (M!, M 2 , M 3 ) - M 5 ) 

+-^= r | (d a M 2 )(d a M 2 ) - (d a M 3 )(d a M 3 ) \ + O f — ) , 
12rM il J \ r J (26) 

<9 t M 6 - ^ (5M 5 - M 8 + M 9 ) + <9 y M 4 J- (M?(M U M 2 ,M Z ) - M 6 ) 
ftM 7 + <9 X M 4 + ^ (5M 5 + M 8 - M 9 ) =± (m^^, M 2 , M 3 ) - M 7 ) , 
<9 t M 8 + d x (M 2 + 3M 6 ) + ^ (M 3 + 3M 7 ) =i (M 8 cq (M!, M 2 , M 3 ) - M 8 ) 

3 ^ J (d a M 2 )(d«M 2 ) + (9jtf 3 )(d a M 3 ) 1 + O ( — ) , 



12 tMi i \ r 

<9 t M 9 - ^ (2M 2 + M 6 ) - ±d y (2M 3 + M 7 ) = l - (m?(M u M 2 , M 3 ) - M 9 ) 



' iA ~ '(,9 a M 2 )(9 a M 2 ) + ( ( 9 a M 3 )( ( 9 a M 3 )l + ^ 



12rMi I \ \ r J 

Thus, we are set up to derive the hydrodynamic equations as the appropriate 
limit of the filtered kinetic system (26). In passing, we note that different mo- 
ments relax with different effective relaxation time scales, because the subgrid 
terms are not the same for all kinetic moments. 



4 Hydrodynamic limit of the filtered kinetic theory 



In the kinetic equation we have a natural length scale set by Knudsen number 
Kn (1). The Navier-Stokes dynamics is obtained in the limit Kn < 1. By 
filtering the kinetic equation we have introduced a new length scale as the 
size of the filter A. The hydrodynamic equations produced by the filtered 
kinetic equation will depend on how A scales with the Knudsen number. In 
order to understand this issue, let us look at the filtered equation for one of 
the moments (26) in the non-dimensional form. In order to do this, let us 
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introduce scaled time and space variables, 

x' = - 

, tc s ( 27 ) 

* = T' 

where c s = l/y/3 for the present model. Let us also specify Knudsen number 
in terms of the relaxation time r: 

where, v is the kinematic viscosity, v = tc\ in the present model. Then, for 
example, the filtered equation for the moment M4 reads: 



d t ,M A + ^d x , (2M 3 + M 7 ) + ^-dy (2M 2 + M 6 ) = ^ (mP(M 1: M 2 , M 3 ) 
+ , o tA^tt 1 (3^0 (d x 'Ms) + (^M 2 ) (^M 3 ) I + O 1 ^ 1 



12KnL 2 Mi v yv y v y yv y \L 4 Kn 



We see that the in absence of the filter (A = 0), the usual situation of a 
singularly perturbed kinetic equation is recovered (and this results in the 
Navier-Stokes equations in the first-order Chapman-Enskog expansion). Let 
us consider the following three possibilities of dependence of A on Kn: 

• If A/L ~ Kn°, then we do not have a singularly perturbed equation in (29) 
anymore. That is, the filter is too wide, and it affects the advection terms 
in the hydrodynamic equations. 

• If A/L ~ Kn, then we do have a singularly perturbed system. However, the 
subgrid terms are of order Kn 2 , and they do not show up in the order Kn 
hydrodynamic equation. In other words, the filter is too narrow so that it 
does not affect hydrodynamic equations at the viscous time scale. 



• Finally, there is only one possibility to set the scaling of filter-width with 
Kn so that the system is singularly perturbed, and the subgrid terms of the 
order A 2 contribute just at the viscous time scale. This situation happens 




Note that, with the scaling (30), all the higher-order terms (of the order 
A 4 and higher) become of the order Kn and higher, so that they do not 
contribute at the viscous time scale. 



Once the scaling of the filter-width (30) is introduced into the filtered moment 
equations (26), the application of the Chapman-Enskog method [10] becomes 
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a routine. We write: 



d t = dj 0) + Kn<9| iJ + 0(Kn 2 ), 



(i) 



(31) 



and for i — 4, . . . , 9: 



Mi = M-\M U M 2 , M 3 ) + KnA/f ) + 0(Kn 2 ). 



(32) 



The hydrodynamics equations at the the order 0(1) are the Euler equations: 



8^ Mi 
d t (0) M 2 = 

9 t (0) M 3 = 



-9 X M 2 - d y M 3 , 
-d x lM lCs + 



M2M3 
" * V Mi 



/ 2 M 3 M 3 



(33) 



Note that no subgrid terms appear at this time scale in the hydrodynamic 
equations (33). This means that large scale motion, even after filtering, is 
dictated just by the conservation laws. Zero-order time derivatives of the non- 
conserved moments are evaluated using the chain rule: 

d^MrCMu^Ms) = ^$f 0) Mi + ^M° } M 2 + ^^ (0) M 3 . (34) 

t , \ 1, 2, i) t m ^ t m ^ t S V ) 



In particular, to the order u : 



di 0) M c 4 \M 1 ,M 2 ,M 3 ) = 0, 
dl 0) MP(M 1 ,M 2 ,M 3 ) = 0, 



M 2 M 2 



a M 2 M 3 



d^Affm, M 2 , Ms) = 9, (M^c 2 . + ^ 

^M-^, M 2 , M 3 ) = + 9, (m iC 1 + ^ 

d^M^M^M^Ms) = 0, 

^ffi.I^Is) = [^M 2 + y M 3 ' 



(35) 



At the next order O(Kn), correction to locally conserved moments is equal to 
zero, 



(i) 



M 



(i) 



M 



(i) 



0, 
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whereas corrections to the non-conserved moments, Mf \ i = 4, . . . , 9, are 
obtained by substituting Eq. (31), and Eq. (32) in Eq. (26) and eliminating 
the zeroth order time derivatives using Eq. (35): 



= -Lc s 



d x M 3 + d y M 2 



+ 



iwf ) = -2 Lc s 



m£ ] = -Lc s 



M { j ] = -Lc s 



d x M 2 - d v M 3 



+ 



12 Mi 
L 2 



d x 



ml) 

Mi 



+ .9, 



6 M 2 M 3 

X = h Oy 

Mi y 



M { * ] = QLc q 



M ( g ] = 6Lc s 



d x M 2 + d y M 3 
d x M 2 + d y M 3 



12 Mi 

QM~ 2 M 3 
Mi 

ml) 

Mi 
3L 2 



(d x M 2 ) (d x M 3 ) + (d y M 2 )(d y M 3 )j, 

d x M 2 f - (d x M 3 ) 2 + (d y M 2 ) 2 - (d y M 3 y \, 



(36) 



12 Mi 

3L 2 
12 Mi 



(d x M 2 ) 2 + (d y M 2 ) 2 + (d x M 3 ) 2 + (d y M 3 ) 2 \, 
(d x M 2 ) 2 + (d y M 2 ) 2 + (d x M 3 ) 2 + (d y M 3 y 



and the first-order time derivative of the conserved moments are: 



a, (1, M, = o, 

3, (1, M 2 = -4 (iM« + 1m»' - ijff ) - 3 s Mi", 

These equations shows that the viscous term and the subgrid term both appear 
as the O(Kn) contribution. We remind that no assumption was made about 
relative magnitude of the subgrid term as compared with the viscous terms. 
The only requirement that is set on the subgrid scale term is that they appear 
at the viscous time scale rather than the time scale of the advection. We can 
write the complete hydrodynamics equation by using Eq. (31), Eq. (32), Eq. 
(33), Eq. (36), and Eq. (37) to obtain the hydrodynamics equations correct up 
to the order 0(Kn 2 ) (at this stage one recovers the Navier-Stokes equations 
using the unfiltered kinetic equation). In the next section, we shall see how 
the subgrid scale terms affect the Navier-Stokes description. 



(37) 
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5 Hydrodynamic equations 



The final set of hydrodynamics equation, valid up to the order 0(Kn 2 ) is: 



d t p + d x pu x + d y puy = 
dtWx~ + d x P xx + d y P xy = 

dtf>Uy + X P X y + dyPyy = 0, 



where 



(38) 



■n . PU X _~ 

P H — 2vpb x 



xy 



yy 



PUx PUy 
P 



- 2vpS, 



xy 



+ P SG 

1 II ) 



~ 1 xy i 



P+^-2UpSyy 



~ 1 vv 1 



(39) 



with P = pc 2 , as the thermodynamic pressure and 



p sg _ KnL2 P 



<? — o 



5*7/3 + ^7/3 



v pL 
12cT 



9 — O 



5*7/3 + ^7/3 



Thus, we have obtained a closed set of hydrodynamics equations, for appro- 
priate choice of filtering width. This set of equations, written in the nondi- 
mensional form up to the order, u 2 is: 



d a u a = 0, 



(40) 



d t (u a ) + dp (u a up) = - d a P + 2Kadp (S a ^ 

~ { - ^7) { s iP + ^7/3) } • 

Note that the pressure appearing in the momentum equation is not the ther- 
modynamic pressure anymore, but needs to be computed from the incompress- 
ibility condition (Eq. 40) [15]. Significance of the Knudsen number appearing 
in the equation (41) is explained below in section 6. 

Similar to the case of the Navier-Stokes equation, the subgrid model (5) enjoys 
the consistent derivation from the kinetic theory. We should remind here again 
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that the result is independent of the particular kinetic model used for the 
derivation. Any kinetic model which recovers the Navier-Stokes equations in 
the hydrodynamic limit will lead to the same result. 



6 Discussion and conclusion 

Now we shall summarise the results obtained in the present work and their 
limitations: 

• It is possible to derive rigorously a coarse-grained closed set of equation for 
hydrodynamics, a long cherished goal in turbulence modeling. 

• The scale-separation present in the kinetic theory provides a natural way 
to obtain coarse models. 

• Arbitrary choice of filter-width is not allowed. 

• In this work, we have shown that the operation of solving the Boltzmann 
equation (Chapman-Enskog expansion) and coarse-graining (filtering) do 
not commute. In the usual procedure of producing filtered hydrodynamic 
equation, the filtering is done on the solution of the Boltzmann equation 
(Navier-Stokes equations), which leads to closure problems. On the other 
hand, reversing the order of these two operations, provides a physical mean- 
ing to the filtering width and produces a closed set of equations in the 
hydrodynamic limit. 

• As the smallest length scale needed to be resolved in the new set of equa- 
tions is A ~ Kn 1 / 2 ~ Re^ 1 / 2 , the cost of numerical computation reduces 
drastically as compared to the fully resolved simulation of the Navier-Stokes 
equations. We can get an estimate of this gain as follows: The smallest scale 
needed to be resolved in the numerical simulation is proportional to the 
Re~ 1/ ' 2 rather than Re~ 3 / 4 (Kolmogorov scale). This changes the scaling of 
number of degrees of freedom in a three-dimensional simulation with Re 
from Re 9//4 to Re 3//2 (number of grid point is oc 5x~ 3 , where 5x is the grid 
spacing). Further, a rough estimate of the scaling of the cost of time inte- 
gration with Re is Re 3//4 (number of time steps is oc 5x~ r ) in the case of 
fully resolved simulation of the Navier-Stokes equations [1]. However, in the 
present case this scaling will be Ma -1 . This happens because any numerical 
scheme has to take a time-step dictated by the sound speed (or an analog of 
it pertinent to the discretization in time chosen). Thus the scaling of the to- 
tal cost of computation with the Re changes from Re 3 for the Navier-Stokes 
equations to Re 3//2 for the present equation. 

• In the above estimation, the Mach number Ma appearing in the equation 
for the filter-width A (A 2 oc Re -1 ) was not taken into the account. This 
is justified as long as Re is sufficiently large and the Mach number Ma is 
not zero. An acceptable limit for the incompressible limit of the Navier- 
Stokes equation is Ma ~ 0.1 (for example, most of the lattice-Boltzmann 
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simulations of the incompressible Navier-Stokes uses Ma ~ 0.05 — 0.1). 
Let us consider the case when the number of grid points in each direction 
is 1024, then a fully resolved simulation of the Navier-Stokes equations is 
possible by taking the Reynolds number as Re ~ O(10 4 ), while a fully 
resolved simulation using the present subgrid model is possible by taking 
the Reynolds number as Re ~ O(10 5 ), for Ma ~ 0.07. 

• One interpretation of the subgrid scale terms is that the removal of the 
small-scales in the kinetic picture appears naturally as the force term. As the 
extra terms appearing in the evolution equations for non-conserved variables 
can also be generated in the unfiltered kinetic equation by an appropriate 
choice of the external force field (dependent on the position as well as molec- 
ular velocity). Thus at least formally, we can find a force-field which will 
act like a filter and remove the small scales of motion present in the kinetic 
equation. Thus the physical meaning of the filtering (a purely mathemati- 
cal operation) at the kinetic level is the application of some self-consistent 
mean-field force which removes the small scale of the motions from the ki- 
netic equation. The technical advantage of the search for a mean-field force 
(appearing in the filtered kinetic equation), rather than an effective viscos- 
ity term (attempts to absorb subgrid scale contributions in the viscous term 
of the Navier-Stokes equations) is that the one does not have to deal with 
the difficult question of what to do with the nonlinearity and nonlocality 
present in the convective term of the Navier-Stokes equations. 

Finally, let us mention some further possible directions of study: 

• From a practical standpoint, a major goal of going beyond a Navier-Stokes- 
based coarse graining, is to make a (filtered) kinetic theory work at possi- 
bly large ratios A/L . Indeed, successful kinetic subgrid models have been 
known empirically for some time [14]. However, it is unclear why kinetic 
subgrid scale models work better than models of the Navier-Stokes equa- 
tions (for example, a recent comment on the kinetic models of turbulence 
was: "Whether the approach can be supported by rigorous theory remains 
to be shown" [16]). Our analysis shows that everything is eventually ruled 
by the smallness of the Knudsen number, a well defined smallness param- 
eter present in the kinetic theory. This is the first rigorous step in the ki- 
netic modeling of the turbulence. For example, the choice of the filter-width 
(2), based on the integer-power (standard) Chapman-Enskog analysis is a 
conservative estimate only. In order to achieve a subgrid model between 
the advection time scale (Kn°) and the viscosity time scale (Kn 1 ), that is, 
A/L ~ Kn 7 with < 7 < 1/2 requires a generalization of the Chapman- 
Enskog method to noninteger series expansion in Knudsen number. This 
interesting possibility needs to be studied separately. The application of 
the method of the invariant manifold [17], which does not require Knudsen 
number Kn to be small, on the filtered kinetic equation for < 7 < 1/2 
is a possible extensions of the present work. The possibility of doing exact 
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Chapman-Enskog expansion [18,19], also need to be investigated further. 
Another possible extensions is the use of "renormalization-group" ideas by 
applying several filters of increasing filter-widths. 

• When the discrete-velocity kinetic theory of section (2) is appropriately 
discretized in time and space, one arrives at the so-called entropic lattice 
Boltzmann method [2,3] (ELBM). In the ELBM, the thermodynamic stabil- 
ity (Boltzmann's if-theorem) is maintained by the discrete-time if -theorem 
[20,21,22,23] which results in unconditionally stable simulation algorithm for 
hydrodynamics. It was argued [21] that ELBM is a built-in subgrid model. It 
would be interesting therefore to establish a closer relation between ELBM 
and the present theory. 

• Filtering the kinetic equations as above can be applied to a wide class of 
kinetic theories with a well-defined separation of time scales enrich existing 
resolved macroscopic models with physically sound subgrid contribution (for 
example, the kinetic equations for the granular flows [24]). 

To conclude, the presented coarse-grained equations are the first rigorously 
derived subgrid model. Effectiveness of the model in practice needs to be 
investigated further numerically. Work in this direction is currently in progress. 
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